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Abstract 

We discuss the connection between the random matrix approach to dis- 
ordered wires and the Calogero-Sutherland models. We show that different 
choices of random matrix ensembles correspond to different classes of CS 
models. In particular, the standard transfer matrix ensembles correspond to 
CS model with sinh-type interaction, constructed according to the C n root 
lattice pattern. By exploiting this relation, and by using some known proper- 
ties of the zonal spherical functions on symmetric spaces we can obtain several 
properties of the Dorokhov-Mello-Pereyra-Kumar equation, which describes 
the evolution of an ensemble of quasi one-dimensional disordered wires of 
increasing length L. These results are in complete agreement with all known 
properties of disordered wires. 
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1 Introduction 



During last years an increasing interest has been attracted by the physics of quan- 
tum electronic transport in disordered wires [Jjj]. One of the main reasons for this 
interest lies in the high degree of universality of some experimental observations. 
In particular it was shown that the variations of the measured conductance as a 
function of the magnetic field or of the Fermi energy are independent of the size 
and degree of disorder of the sample and have a variance always of order e 2 /h. This 
phenomenon is usually known as universal conductance fluctuations (UCF). Its uni- 
versality suggests that disordered wires could be described by some relatively simple 
Hamiltonian, independent of the particular model or disorder realization. This ap- 
proach was pioneered by Imry and developed by Muttalib, Pichard and Stone ||, 
who suggested to describe UCF by constructing a Random Matrix Theory (RMT) 
of quantum transport, in analogy to the Wigner-Dyson RMT for nuclear energy 
levels. However it was soon realized that, besides the obvious analogies, there are 
some relevant differences between the Wigner-Dyson ensembles (WDE) of random 
matrices and those constructed to describe the physics of disordered wires (which 
we shall denote in the following as Transfer Matrix ensembles (TME)). The most 
interesting feature of these TME is the presence of a Fokker-Plank type evolution 
equation, known as DMPK equation (see below) which can be considered the equiv- 
alent in the context of TME of the Brownian motion approach to WDE suggested 
by Dyson (for a discussion of these analogies see M). 

The aim of this contribution is to show that the relationship between WDE and 
TME (and the parallel one between Brownian motion and DMPK equation) can 
be well understood by exploiting their equivalence with the so called "Calogero- 
Sutherland" (CS) models || which are quantum integrable models describing a set 
of N particles moving along a line (see below). In particular it will be shown that 
the Schrodinger equation of the CS models is equivalent to the DMPK equation and 
that WDE's, TME's and the S-matrix ensembles described in [[J correspond to dif- 
ferent realizations of the CS models. We shall show that the common, underlying, 
mathematical structure of all these models is the theory of Laplace-Beltrami oper- 
ators on Symmetric Spaces. By using some recent results on the eigenfunctions of 
these operators, which are known as "zonal spherical functions" we shall construct 
exact (for the unitary ensemble) or asymptotic (in the other cases) solutions of the 
DMPK equation and describe several physical properties of disordered wires. 

This contribution is organized as follows: after a short introduction to disorder 
wires and to their TME description (sect. 2) we shall discuss the DMPK equation 
(sect. 3). In Sect. 4 we shall give a short introduction to the CS models and show 
the anticipated equivalence with the TME. Sect. 5 will be devoted to the solution 
of the DMPK equation and sect. 6 to some concluding remarks. 
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2 Disordered wires 



The peculiar feature of disordered wires is that in these systems, at low enough 
temperature, the phase coherence of the electron's wavefunctions is kept over dis- 
tances much larger than the typical mean free path, thus allowing to study several 
non-trivial quantum effects. These phenomenon can be studied only at low tem- 
peratures where the inelastic electron-phonon scattering (which changes the phases 
of electrons in a random way) becomes negligible and resistivity is completely dom- 
inated by the scattering against random impurities, which is elastic and changes 
the phases in a reproducible, deterministic way. The natural theoretical framework 
to describe these systems is the Landauer theory |7j which assumes the electrons 
in thermal equilibrium with the various chemical potentials in the leads. The dis- 
ordered wire, with its impurities, is then regarded as a scattering center for the 
electrons originating from the current leads and the conductance G is proportional 
to the transmission coefficients of the scattering problem. Within this approach 
Fisher and Lee || proposed the following expression for the conductance in a two- 
probe geometry (namely a finite disordered section of length L and transverse width 
W, to which current is supplied by two semi- infinite ordered leads): 

G = G Ti(tt^=G Y,T n , G = ^f. (1) 

where t is the N x N transmission matrix of the conductor and 7\ , T 2 • • • are the 
eigenvalues of the product ttf and are usually called transmission eigenvalues. N is 
the number of scattering channels at the Fermi level. N depends on the width of 
the wire and even in the narrowest metal wires it is of the order of N ~ 10 4 — 10 5 
so that for metal wires a large N approximation will give in general very good 
results (notice however that semiconductor microstructure with very low values of 
N can be constructed and studied). In the following we shall often refer to the 
dimensionless conductance g, defined as g = G/Gq. The transmission matrix t is a 
component of the 2N x 2N scattering matrix S which relates the incoming flux to 
the outgoing flux: 

s = ( 1 1> ) • < 3 > 

where I,0,I',0' are iV component vectors which describe the incoming and outgoing 
wave amplitudes on the left and right respectively, and r is the N x N reflection 
matrix. Current conservation implies 

|/| 2 + |/'| 2 = \0\ 2 + \0\ 2 , (4) 

which is equivalent to the requirement of unitarity: S G U(2N). However it turns 
out that for the problem that we are studying a much better parametrization is 
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given by the transfer matrix M which is defined as: 





(5) 

The relation between the transfer matrix and the transmission eigenvalues be- 
comes clear if one constructs an auxiliary matrix Q defined in terms of M as follows: 

Q = - [M^M + (M^M)- 1 - 2] (6) 

The eigenvalues {Aj} of Q are non-negative and can be related to the transmission 
eigenvalues Tj by 

Xi = {l-T^/Ti (7) 

The physics of the disordered wires that we are studying will be completely de- 
scribed if we can obtain the probability distribution P({Aj}) for the eigenvalues Aj 
(and consequently for the Tj's). 

To this end let us first study some general properties of the transfer matrix 
which are direct consequences of the physical symmetries of the problem. First 
of all, it is easy to see that the same flux conservation constraint eq.(f|) discussed 
above implies in this case the conservation of a hyperbolic norm: 

M ] T, Z M = S z (8) 

with 

(9) 

where and 1 are the zero and unit N x N matrices. As a consequence of (|8|), 
M e SU(N, N). The ensemble of transfer matrices defined in this way is usually 
called "unitary ensemble" (ensemble Ha in the notation of ref. pi). 

If the system is also invariant under time reversal symmetry, M must satisfy a 
further constraint: 

.1/ •>;,..!/ = E x (10) 

with 

/ n i \ 

(11) 

It is possible to show that the joint application of @ and ( |10|) implies M 6 
SP(2N, R) ||10|| . The ensemble of transfer matrices defined in this way is usually 





called "orthogonal ensemble" (ensemble I in the notation of ref. 

From an experimental point of view it is very simple to control the time reversal 
symmetry which is eliminated by the application of an external magnetic field. 

If the disordered wire contains "magnetic impurities" , namely if the spin-orbit 
interaction in the scattering against impurities becomes important then the spin- 
rotation symmetry (which was implicitly assumed in all the above discussion) is not 
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any more conserved. In this case the two spin components of the electrons must be 
treated separately. Each one of the input and output vectors J, I', O, O' becomes a 
collection of iV spinors (one for each channel) and each spinor contains the two spin 
degrees of freedom. Hence M is in this case a 4N x 4iV complex matrix. If only flux 
conservation is imposed (time reversal symmetry broken) then M G {/ (2N, 2N) and 
we find again the unitary ensemble described above with the only change: N — ► 2N 
(this difference was kept explicit in ref. |J, where this case was denoted as lib 
to distinguish it from the spin- rotation symmetric unitary ensemble Ha). If time 
reversal symmetry is conserved (namely if there are magnetic impurities, but no 



external magnetic field) then one must impose the following constraint on M [p 
(analogous of that of eq.flToD) 



with 



v A 

where A is a 2N x 2N block diagonal matrix 



M* = KMK T (12) 



K= l X A • (13) 




A = <7l, „=[ " ' . (14) 



This constraint implies that M G SO*(4N) and defines the so called symplectic 
ensemble (ensemble III in the language of ref |§). 

The second step in order to construct the probability distribution for the Aj's 
comes from the identification of the Aj's themselves as the relevant physical de- 
grees of freedom of the system. This identification has some very interesting group 
theoretical consequences. In fact the choice of the Aj's as relevant physical param- 
eters induces the following parametrization for M for the orthogonal and unitary 
cases [HOl: 













,(3) 












M = ('" 

X ( "o u« J S UTV (15) 

where A is a iV x real, diagonal, matrix with entries the eigenvalues Ai, A2, ■ • ■ A^ 
in both cases. The u^ l \ (i = 1, 2, 3, 4) are 4 independent N x N unitary matrices in 
the unitary case while in the orthogonal case they are constrained by the relations: 

u m = u w^ u w = u m* (i 6 ) 

In the symplectic case we have again the same parametrization, if the various 
matrices are written in terms of quaternions. Thus in this case A is a. N x N 
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quaternion real, diagonal matrix and the u'^'s are N x N quaternion, unitary 
matrices which again obey the constraint (fL6|). 

In this parametrization we recognize a G/H coset structure, where G is the 
group to which the transfer matrix belongs and H is the group to which the U and 
V matrices belong. The coset structure is immediately evident if we notice that any 
transformation M —>■ M' = WMW~ l , with W e H gives again a transfer matrix 
which can be decomposed as M' = U'TV with U' = WU, V = VW~ l and the 
same matrix V. So the physically relevant parameters \ are left unchanged by such 
transformation and thus belong to the coset space G/H. These cosets are listed in 
Tab.l for the three ensembles in which we are interested. 



G 


H 


T 


S-R 


Sp(2N, R) 


U(N) 


y 


y 


SU{N, N) 


SU(N) ® SU(N) ® U(l) 


n 


y 


SO* {AN) 


U(2N) 


y 


n 



Table 1: G/H cosets for the orthogonal (first line), unitary (second line) and 
symplectic (last line) ensembles. In the first two columns the group G and the 
subgroup H. In the last two columns the status of the time reversal (T) and spin 
rotation (S-R) symmetries respectively (y= conserved, n=broken). 

However this is not the end of the story. Looking at the three particular real- 
izations of the pair G, H listed in tab.l we see that in all the three cases the cosets 
are actually symmetric spaces (see tab. 2). What is more important, we recognize in 
the parametrization (|l"5|) the so called "spherical decomposition" of those symmetric 
spaces (see for instance This tells us that the Aj's play the role of (general- 

ized) radial coordinates in G/H and implies that the Aj's are also invariant under 
(generalized) angular transformations in G/H and that only the radial projection 
of any given dynamical operator will influence their dynamics. 

At this point the only remaining step is to impose some dynamical principle so 
as to obtain an "equation of motion" for the probability distribution of the Aj's. 



3 The DMPK equation. 

This program was completed during the eighties, at least in the case of quasi one- 



dimensional wires, by Dorokhov ||13|| , and independently by Mello, Pereyra, and 
Kumar [|l(J (for j3 — 1) by looking at the infinitesimal transfer matrix describ- 
ing the addition of a thin slice to the wire. The resulting evolution equation for 
the eigenvalue distribution P({Aj}, s) is usually known as Dorokhov-Mello-Pereyra- 
Kumar (DMPK) equation. The only assumptions which are needed to obtain this 
equation are first that the conductor must be weakly disordered so that the scat- 
tering in the thin slice can be treated by using perturbation theory and second 
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that the flux incident in one scattering channel is, on average, equally distributed 
among all outgoing channels. It is exactly this second assumption which restricts 
the DMPK equation to the quasi 1-d regime, where finite time scale for transverse 



diffusion can be neglected. The results of [10] were then generalized to (3 = 2,4 in 



Refs. [0, |TJ. The DMPK equation is: 

dP 

^ = DP, (17) 

where s is the length L measured in units of the mean free path I: s = L/l and D 
can be written in terms of the Aj's as follows: 

D = - £ ^A,(l + A,) J{\)£- J(A)-\ (18) 

with 7 = (3N + 2-/3. (3 G {1,2,4} is the symmetry index of the ensemble 
of scattering matrices: (3 = 1 for the orthogonal ensemble, f3 = 2 for the unitary 
ensemble and (3 — 4 for the symplectic one, in full analogy with the well know 
Wigner-Dyson classification. J({A n }) denotes the Jacobian from the matrix to the 
eigenvalue space: 

= n '.\ x j - a /- ( i9 ) 

i<j 

There is however a completely independent, and very elegant, way to obtain 
the DMPK equation. Let us assume as dynamical principle to obtain an equation 
of motion for P(X) the simplest possible choice compatible with the constraints 
described in sect. 2. That is, let us assume that, as L (the length of the wire) 
increases, the matrix M freely diffuses in the G/H space going from the perfectly 
conducting limit (L = 0) (which plays the role of initial condition for this evolution 
equation) to the insulating, localized limit (L = oo). In general a free diffusion 
in G/H is described by the Laplace-Beltrami operator of G/H. However, since 
the A's are the radial coordinates of G/H, their behaviour as a function of L will 
only depend on the radial part B of the Laplace-Beltrami operator. The resulting 
equation is: 

OP 

— = a BP, (20) 

where a is a (for the moment undetermined) diffusion constant and B is defined as 
follows: 

B =^±£-mr£- t , (2D 

where we have chosen the following parametrization for the radial coordinates of 
the manifold: Aj = sinh 2 Xj, and 

£(x) = Y[ I sinh 2 Xj — sinh 2 ^-i I ^ TT I sinh2xj| 5 . (22) 

i<3 i 
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It is now only matter of straightforward algebra to recognize that B and the 
DMPK operator D are related by: 



D 



- [£(*)] 2 B Ux)}-\ 
27 



(23) 



thus allowing, through a suitable choice of a and normalization of P(X) to iden- 
tify eq. (fL7D and eq.(^0|). This identification was first recognized by Hiiffmann ||15| , 
and has been recently discussed in [[16| and fTTj] . 

Eq.(^) can be considered the analogous, in the context of the TME's (hence 
for symmetric spaces of negative curvature) of the Brownian motion approach in 
the case of S-matrix ensembles (which as we shall show below are characterized, by 
the same symmetric spaces, but with positive curvature). 

An important and unexpected property of the DMPK equation is that if 3 = 2 
the various Aj can be decoupled. This was recently realized by Beenakker and 



Rejaei |18] who showed that the DMPK equation can be rewritten as Schrodinger- 
like equation (in imaginary time) for a set of A" interacting fermions. The mapping 
was obtained by setting: A„ = sinh 2 x n , and by making the following substitution 



P({Xn},s) = 

In this way the DMPK equation becomes exactly equivalent to: 



(24) 




27 



i<j 



dxf sinh 2 2x, 



sinh 2 2xj + sinh 2 2xi 
(cosh 2xj — cosh 2xi) 2 



N 3 3 2 

U = A^(A^ - 1)- - N(N - 1)(N - 2) — . 

27 7 67 



(25) 



(26) 



(27) 



By choosing 3 = 2, the remaining interaction terms among the x n disappear, the 
equation can be decoupled and can be solved exactly 



This equivalence with a Schrodinger equation is another feature of the DMPK 
equation which has a natural explanation in the framework of the Calogero- Sutherland 
models. 



4 Calogero- Sutherland Models 

These models describe A^ particles on a line, identified by their coordinates {xi}, i = 
1 • ■ • N, interacting (at least in the simplest version of the models) with a pairwise 
potential f(xi,Xj). Several realizations of this potential have been studied in the 
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literature (for a comprehensive review see ref. |JI9|), but in the following we shall 
mainly be interested in only two realizations. The sin-type CS model, in which 
f(xi,Xj) = 1/ sin 2 (xj — xj) and the sinh-type for which f(xi,Xj) = 1/ sinh 2 (xj — Xj). 

The most relevant feature of these models is that (under particular conditions 
discussed below, see eq.(^8|) and ([34])) they have N commuting integrals of motions, 
they are completely integrable and their Hamiltonian can be mapped into the radial 
part of a Laplace-Beltrami operator on a suitable symmetric space. In particular we 
have spaces with negative curvature for the sinh models and of positive curvature 
for the sin-type ones. 

In the original formulation of the CS model, the interaction among the particles 
was simply pairwise 0. But it was later realized that the complete integrability 
of the model had a deep group theoretical explanation, that the simple pairwise 
interaction was the signature of an underlying structure: namely the root lattice 
of the Lie algebras A^, and that all the relevant properties (complete integrability, 
mapping to a Laplace-Beltrami operator of a suitable symmetric space) still hold 
for potentials constructed by means of any root lattice canonically associated to a 



simple Lie algebra [19]. Let us see more precisely how this construction works. 

Let us call V the N dimensional space defined by the coordinates {xj} and 
x = (xi, • • -xn) a vector in V. Let R = {a} be a root system in V, and R + the 
subsystem of positive roots of R. Let us denote with x a the scalar product (a,x). 
Then the general form of the CS Hamiltonian is 

1 N d 2 g 2 
^ = -2§^ + a g + ^?fc) (28) 

where the couplings g 2 are the same for equivalent roots, namely for those roots 
which are connected with each other by transformations of the Coxeter group W of 
the root system. To clarify this rather abstract definition let us see two examples, 
obtained using the root lattices A N and Cn (in the following {e^, • • ■ e^} denote a 
canonical basis in the space R n ). 

An'- This root system is obtained by taking a hyperplane in R^" 1 " 1 for which 
X1 + X2 + ■ ■ ■ xn+i = 1. Then the root system R is given by: R = {ej — ej, i 7^ 
j}. In this case W is the permutation group of the set {ej}. The corresponding 
Hamiltonian is: 

1 N d 2 g 2 
^^g^ + gsinh 2 ^-^) (29) 

This is the model originally considered in || 

Cn- This root system is R = {±2ej, ± ej ± ej, i 7^ j}, in this case W is the 
product of the permutation group and the group of transformations which 
change the sign of the vectors {ej}. The corresponding Hamiltonian is: 
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N 



° ^ ;)2 - Ysinh 2 (2x i 



2 ^ <9 2 X,; 

+ 



4? 



2 2 



+ ^ \sinh 2 (xj — x,-) ' sinh 2 (a; i +Xj)/ 

This is the model which we shall study in the following. 

By using simple identities among hyperbolic functions eq. 
as follows: 



H 



+ c 



2 ^ sinh 2 (2a;j) + sinh 2 (2x :) ^ 



i<j 



(cosh(2xj) — cosh(2;Cj))' 



with c an irrelevant constant. By setting 



(30) 



can be rewritten 



(31) 



(32) 



we see that eq. fl3~0f) coincides (apart from the overall factor I/7) with TC in eq. (p6|) . 

As we mentioned above, the relevant feature of the CS hamiltonian (^) is that 
it can be mapped into the radial part B of a Laplace-Beltrami operator (see for 
instance Appendix D of ref. |19[]) of a suitable symmetric space 



H = £(x) 



(33) 



with B defined by eq.(|2~T|) and p a constant term which we shall neglect in the 
following. The particular symmetric space is uniquely fixed by the root lattice 
underlying the CS Hamiltonian and by the coupling constant g a In fact it is well 
known that all the irreducible symmetric spaces of classical type can be classified 
with essentially the same techniques used for the Lie algebras. They fall into 11 
classes labelled by the type of root system and by the multiplicities of the various 



roots Jl2|. Some of these spaces (those relevant for our discussion) are listed in 
tab. 2 and 3 with their root multiplicities. 



These multiplicities are related to the coupling constants by [19 



2 m a {m a -2) 2 
9 a = o a 



(34) 



where \a\ is the length of the root a and m a its multiplicity. Only for these special 
values of g 2 the mapping (^) is possible. For the three ensembles in which we are 
interested we have m a = f3 for the short roots (those of the type {±ej ± e,-}) and 
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G 


H 


P 


SL(N, R) 


SO(N) 


1 


SL{N, C) 


SU{N) 


2 


SU*(N) 


USp{2N) 


4 



Table 2: Irreducible symmetric spaces of type An- In the first two columns the 
group G and the maximal subgroup H which define the symmetric space. All these 
spaces are labelled by the Dynkin diagram An-\- In the last column the multiplicity 
P of the roots 



G 


H 


P 


7? 


Sp(2N, R) 


U(N) 


1 


1 


SU(N, N) 


SU(N) g> SU(N) ® U(l) 


2 


1 


SO*{4N) 


U{2N) 


4 


1 


USp{2N, 2N) 


USp{2N) ® USp{2N) 


4 


3 


Sp{2N, C) 


USp(2N) 


2 


2 



Table 3: Irreducible symmetric spaces of type Cn- In the first two columns the 
group G and the maximal subgroup H which define the symmetric space. All these 
spaces are labelled by the Dynkin diagram Cn- In the third column the multiplicity 
P of the ordinary roots ot Cn- In the last column the multiplicity rj of the long 
root. 

m a = 1 for the long roots (those of the type {2e,}), which if inserted in eq. (|34D 
exactly give the values of eq.(^). In this way we see that the index P has a deep 
group theoretical meaning, since it denotes the multiplicity of the (ordinary) roots 
of the symmetric space in which the transfer matrix diffuses as L increases. The 
mapping described by eq.(|33|) is exactly the one found in ref. 

All the symmetric spaces listed in tab. 2 and 3 are spaces of negative curvature. 
For each one of them there is a counterpart of positive curvature, with all other 
properties (in particular the root lattice structure) unchanged. For instance we 

haV6 ' SL(N,R) SU(N) 

SO{N) ~* SO{N) 1 ' 

In the context of CS models one moves from negative to positive curvature sym- 
metric spaces by changing the sinh-type interaction into the sin-type one. In the 
framework of RMT for disordered wires one has the same change moving from the 
TME to the S-matrix ensembles. In fact, as we mentioned above, the flux conser- 
vation constraint implies that the S matrix belongs to a compact group. At this 
point, depending on the problem in which one is intersted, and consequently, on 
the parametrization which one choses for the (possibly generalized) eigenvalues one 
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can find WDE's, which correspond to sin CS models of type, or the S- matrix 
ensembles described in || which correspond to sin CS model of CV type. 



5 solution of the DMPK equation 

The most important application of the above described equivalence between TME 
and CS models is that, by using some recent results on the CS models one can obtain 
several important properties of the DMPK equation: solve it exactly in the case (3 = 
2 and find approximate asymptotic solutions for (5 — 1,4. This result was obtained 
in |16| and we shall describe here the main steps of that solution. According to 
eq.(p3|) if $fc(x), x = {x%, ■ ■ ■ ,x^}, k = {hi, • • • , kj^} is an eigenfunction of B with 
eigenvalue k 2 , then t;(x) 2 &k{%) will be an eigenfunction of the DMPK operator 
with eigenvalue k 2 /(2 / y). These eigenfunctions of the B operator are known in 
the literature as "zonal spherical functions". In the following we shall use three 
important properties of these functions (see pOf). 

1] By means of the zonal spherical functions one can define the analog of the Fourier 
transform on symmetric spaces: 



/(*) 



Kk)$ k {x) 



dk 



|c(*)| s 



(36) 



(where we have neglected an irrelevant multiplicative constant) and in the three 
cases which are of interest for us: 



|c(A;)| 2 = |A(fc)| 2 n 



r(^) 



with 



|A(*)| S 



n 



r(< 



r(i + 
r(i 



§_ I • fcm fcj 
2 "'" L 2 



f3_ _|_ ■ k m +. 
2 ~T 1 2 



(37) 



(3* 



where T denotes the Euler gamma function. 

2] for large values of x, &k(x) has the following asymptotic behaviour: 

1 



fa) 



c{rk)e i(rk ' x) 



(39) 



where rk is the vector obtained acting with r G W on k. The important feature of 
eq. fl39D is that it is valid for all values of k. 



3] in the case (3 = 2 the explicit form of is known JT9|, |21 

det [gi] 

rii<j l(k 2 — kj) (sinh 2 Xi — sinh 2 Xj) 



(40) 
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where the matrix elements of Q are: 

Ql = F Q(l + ik m ), - ik m ), 1; - sinh 2 % ) (41) 

and .F(a, 6, c; z) is the hypergeometric function. 

Eq.s( |23| , |36| - p8D allow to write the s-evolution of P({x n },s) from given initial 
conditions (described by the function fo(k)) as follows: 

/fc 2 Hk 
fo(k)e-^ h {x)j^. (42) 

By inserting the explicit expression of |c(fc)| 2 and by using the identity: 




we end up with the following general expression for P({x n }, s) with ballistic initial 
conditions (which, due to the normalization of $fc(x), simply amount to choosing 
fo{k) = const): 



P({x n },s) = [ax)} 2 Jdke~^ 



|A(*)| S 



n^tanh(^) (44) 

3 

This expression is rather abstract, but it can be made more explicit by using the 
properties [2] and [3] listed above. In the (3 = 2 case we can insert the explicit 
expression for $fc(x), given in eq.s( [4T)| , f4"l"D , into eq.(fPJ). By using the identity 

P„(z) =F(-u,u+ 1,1; (l-z)/2) (45) 

we exactly obtain (as expected) the solution, found by Beenakker and Rejaei in the 



same case [[18 

In the other two cases (5 = 1, 4, if x is large (and in our framework this means 
x 2 > (2s) /j) we may insert the asymptotic expansion (p9[) into eq.(^). The re- 
sulting behaviour of P({x n },s) will depend on the chosen (metallic or insulating) 
regime for k. Let us look at the two cases separately. 

Insulating regime (k <^1). 
In the k — > limit the T functions in eq. (|39|) can be approximated according to: 

r(@- + iv) 

T{iy) ~ y ^° iy ' ^ G { 1 ' 2 ' 4 }' ( 46 ) 
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Then, by rewriting both the product l\i<j(k 2 — k 2 ), and the sum over the exponen- 
tials in fl39|) as determinants, the integration over k becomes straightforward and 
gives: 



P({x n }, s) = Y\ sinh 2 £j — sinh 2 a;, 



X : 



i<j 



X 



JJ [exp(-s 2 7/(2s))x i (sinh2x i ) 1/2 



(47) 



Ordering the x ra 's from small to large and using the fact that in this regime 
1 < ii <C x 2 <C ■ ■ • < % we can approximate the eigenvalue distribution as 
follows: 



A' 



P({x n },s) = ]Jexp -(7/(2s))(z« - x^ 2 



(4£ 



where x n = ~(1 + j3(n — 1)), in agreement with the result obtained by Pichard,f22 
by directly solving the DMPK equation in this regime. 

Metallic regime (k 3> 1). 
In this case one must use the asymptotic expansion: 



\y\ 2 e 4 



Pe {1,2,4} 



(49) 



The integration over k is less simple in this case and, (to be consistent with the 
regime of validity of eq. (|39|) ) in the resulting expression only the highest powers of 



( x 



must be taken into account. We find: 



P({x n }, s) = Y[ |sinh 2 Xj — sinh 2 x 

i<j 

n [ ex p 



2 _ 2 

Xj x i 



X 



x 2 7/(2s))(xi sinh 2x;) 1/2 



(50) 



In agreement with the exact result of [[18] for (3 = 2 and with the (3 dependence 
found by Chalker and Macedo through a direct integration of the DMPK equa- 
tion. 

According to ref. [B5| and ref. [20 1 , the regime of validity of eq.s(|47|) and (^0) is 
x 2 > (2s)/ j, which gives in in the large N limit x 2 > (2s)/((3N). Notice however 
that eq. (|39"D is only the first term of a series which converges absolutely to $fc(x) 



for all values of k (see sect. 8 of ref. 19 and ref. 201). The coefficients of this series 



can be constructed recursively, thus allowing to study the behaviour of eq. (f44"[) even 
for values of x smaller than the above mentioned threshold. 
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6 Conclusions and Perspectives 



We have shown that there is a deep connection between CS models and the RMT 
approach to disordered wires. In particular we have shown that sinh CS models 
of Cat type correspond to TME's, that sin CS models of Cn type correspond to 
the S-matrix ensembles in the parametrization of ref || and that ordinary WDE's 
correspond to sin CS models of An type. Several properties of the TME's can be 
obtained by exploiting this connection. In particular exact or approximate solution 
of the DMPK equation can be obtained. 

Another interesting application of this correspondence can be found in the in- 
sulating regime. Looking at tab. 2 and tab. 3 we see that the most relevant feature 
of the Cn symmetric spaces with respect of the A N ones is that they are described 
by two critical indices instead of one. This is due to the fact that the Cn Dynkin 
diagrams have two types of roots, each with its particular multiplicity. We shall 
call this second index r\ in the following. This peculiar feature of these ensembles 
is usually ignored because in the weakly localized regime it is only the value of /3 
which matters, and also because the three cases which have been studied up to now 
(the first three lines of tab.l) have the same value of the second index 77 = 1, thus 
leading to the same DMPK equation (with no explicit rj dependence). 

However in the insulating regime also the index rj becomes important and can 



be directly measured by looking for instance at the ratio |24 

var (log g) 2 



(51) 



< log# > r) 

In fact, while in the weakly localized regime all the relevant physical properties 
are completely determined by the level statistics (namely the value of /3), when 
extending the RMT approach to the insulating regime, due to the fact that in 
this case the conductance is dominated by the lowest eigenvalue all the details of 
the chosen RMT model (number of degrees of freedom, possibly the presence of 



new critical indices) become important |24|. This observation could open a new 



interesting field of application of the Calogero- Sutherland techniques in the physics 
of disordered wires 
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